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Abstract 

We report preliminary measurements of the branching fractions for B + — > J/ip(j)K + , B° — ► 
J/ip(j)Kg, B° — > J/ip4>, B° — > J/iprj and 5° — > J/iprj' using 56 million BB events collected 
at the T(4S) resonance with the BABAR detector at PEP-II. We measure branching fractions of 
B{B+ -» J/^(j)K + )= (4.4 ± lA(stat) ± 0.7{syst))x 10" 5 and B(B° -> J/^ifg)=(5.1 ± 1.9(stat) ± 
0.9(syst)) x 10~ 5 , and set upper limits at 90% C.L. for branching fractions B(B — ► J/ip(j))< 
0.95 x 10~ 5 , B(B° -» J/iprj)< 2 -? x !0~ 5 , and ^(B -> J/ipr)')< 6.4 x 10" 5 . 
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1 Introduction 



The Cabibbo-favored transition b — > ccs is well established by observation Q of B decays to a 
charmonium state and a kaon, such as B — > J/ipK and J/tpK* . Recent observations of the decays 
-B — » J/?/;7r |Q] and @ are evidence for the Cabibbo-suppressed transition 6 — ► ccd. The quark 

diagrams for these color-suppressed decays are shown in Figures [l] (a) and (b). We search for B 
meson decays into other final states with charmonium. Since B — > J/^tt is observed, the Cabibbo 
suppressed modes J/^rj and J/ipr]' may exist at a comparable level. A further test is to search for 
quark combinations such as bq — > ccsssq, where the ss quark pairs are produced from sea quarks 
or are connected via external gluons as shown in Figures |2| (a) and (b) . This would be exemplified 
in modes such as B — > J /ip(j)K. The mode J/4>4> is a pure rescattering process, the measurement of 
which can help to resolve the discrete ambiguity in the cos{2[3) measurement with B — > J /ipK* Q. 
In this paper we report a search for B decays into J/ip<ft, J/iprji J/^Pv' \ J/tp4>K + , and J/ip(j)Kg and 
present their branching fractions or upper limits. 

Using a factorization approximation with heavy quarks, A. Deandrea et al. [|| have predicted 
the branching fraction of B — > J/ipr] to be a factor of 3.7 smaller than B° — > J/tpn , corresponding 
to B(B° — > J/ipr]) = 0.54 x 10~ 5 Hj. The L3 Collaboration [0] searched for this mode, found no 
events and set an upper limit, B(B° — > J/iprj) < 1.2 x 10 -3 at 90% confidence level. The mode 
B — > J/ij}(j)K has been observed by the CLEO Collaboration with 10 events in 9.6 x 10 6 -B-B 
pairs with the result — > J/ip^K) = db 1.3^ x 10 -5 . In addition to yielding cc bound 

states, the decay i? — » X(cc) + -RT may provide hybrid charmonium (cc + ^/ne) and the hybrid 
state may ultimately decay into J /ip<f> in the final state J /ip<j)K. No published results exist for the 
modes B — > J/^ and J/tprj'. 



(a) 



(b) 




Figure 1: Quark diagrams for (a) -B — > J/tpK and J/ipK* and (b) -B — ► J/i/jtt and J/ipp. 



2 BABAR Detector and Dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e + e~ storage ring. The complete detector is described in detail elsewhere Q. We briefly describe the 
relevant detector subsystems for the physics analysis in this paper. The BABAR detector contains 
a five-layer silicon vertex tracker (SVT) and a forty-layer drift chamber (DCH) in a 1.5-Tesla 
solenoidal magnetic field. These devices detect charged particles and measure their momentum and 
energy loss. The transverse momentum resolution is cr pt /p t = (0.13 ± 0.01) x p t % + (0.45 ± 0.03)%, 
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q « q 



Figure 2: Quark diagrams for B — > J ~/ip(f)K via (a) strange sea quarks and (b) gluon coupling. 



where pt is measured in GeV/c. Photons and neutral hadrons are detected in a CsI(Tl) crystal 
electromagnetic calorimeter (EMC). The EMC detects photons with energies as low as 20 MeV 
and identifies electrons by their large energy deposit. The EMC energy resolution for photons and 
electrons is a(E)/E = 2.3%/£(GeV) 1 / 4 + 1.9%. The charged particle identification (PID) combines 
SVT and DCH track energy loss measurements and particle velocity measurements by an internally 
reflecting ring-imaging Cherenkov detector (DIRC) of quartz bars circumjacent to the DCH. The 
slotted steel flux return is instrumented with 18-19 layers of planar resistive plate chambers (IFR). 
The IFR identifies penetrating muons and neutral hadrons. 

The data used in these analyses were collected in two periods, October 1999 to October 2000 
and February 2001 to December 2001. The data correspond to a total integrated luminosity of 
~ 51 fb" 1 taken on the T(45) resonance and 6.3 fb _1 taken off-resonance at an energy 0.04 GeV 
below the T(4S) center of mass energy and below the threshold for BB production. This data set 
contains ~ 56 million BB events (N B -g). 

3 Physics Analysis 
3.1 Particle Selection 

This analysis begins with selection of charged particles and photons. All charged particle track 
candidates must have at least 12 DCH hits and pt > 100 MeV/c. The track candidates not associ- 
ated with a K® decay must also extrapolate to a nominal interaction point within \J x 1 + y 2 <1.5 
cm and \z\ <3 cm where the origin is at the interaction point, the z axis is along the electron 
beam direction, the y axis is vertically up, and the x axis points away from the collider center. 
The muon, electron, and kaon candidates must have a polar angle in radians of 0.3 < 0^ < 2.7, 
0.41 < 8 e < 2.409, and 0.45 < 9k < 2.5, respectively. In addition, all charged kaon candidates used 
in this analysis are required to a lab momentum greater than 250 MeV/c. These restrictions keep 
the tracks in regions that are well understood by the PID systems. 
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Photons candidates are identified as hits in contiguous EMC crystals that are summed together 
to form shower clusters and have a minimum 30 MeV shower energy and satisfy certain shower shape 
criteria expected for electromagnetic showers. The variables that describe the shower shape include 



the lateral energy |l0| (LAT) that determines the radial energy profile, and Zernike moment [11| 
(A42) that measures the asymmetry of the cluster shape about its maximum. For electron showers 
the LAT peaks near 0.25 and A42 peaks near zero. All the photon candidates are required to have 
LAT< 0.8. 

Electron candidates are required to have a good match between the expected and measured 
energy loss (dE/dx) and between the expected and measured DIRC Cherenkov angle (9c)- Also 
the measurements of the ratio of EMC shower energy over DCH momentum (E/p), and the number 
of EMC crystals associated with the track candidate must be appropriate for an electron. We define 
very tight (VTE) and loose (LE) electron selection criteria that have efficiencies of 88% and 97%, 
respectively. 

Muon candidates are required to have measurements of several variables that help distinguish 
muons from other charged particles. These measurements are: the EMC energy, the number of 
hit layers in the IFR, the penetration depth expressed in units of interaction length along the 
track's path in the IFR and EMC, the difference between the expected and measured number of 
interaction lengths, the average number of hits per IFR layer, the variance of the distribution of the 
number of hits on each IFR layer, the fraction of hit layers between the innermost and outermost 
layer, the chi-square match of hits in the IFR, and the chi-square match between the IFR and the 
extrapolated DCH track. We combine these variables to form different selection criteria applicable 
in different modes. The criteria are called tight (TM, efficiency 70%), loose (LM, efficiency 86%), 
and very loose (VLM, efficiency 92%). 

Charged kaon candidates are selected based on dE/dx information from the SVT and DCH 
and 9c- A likelihood function that combines all the information is constructed for the kaon, the 
proton and the pion hypotheses. A likelihood ratio test determines if the candidate track satisfies 
the loose kaon selection (LK), very tight kaon selection (VTK) or the not-a-pion selection (NP). 
The SVT, DCH and DIRC information and the likelihoods are used in certain selected momentum 
ranges. The loose and very tight selections have typical efficiencies from 70 to 90%, whereas the 
extremely loose selection, not-a-pion, has > 90% efficiency. 



3.2 Event Selection 

The estimation of the signal and the background employs two kinematic variables; the energy 
difference AE, which is the energy of the B candidate in the T(AS) frame minus the energy of 

the beam particle, E^^ m , and the energy substituted mass Mes which is J (-E^m) ~~ (Pb M ) 2 1 
where Pjj is the momentum of the B candidate in the T(4S) frame. Typically these two weakly 
correlated variables form a two dimensional Gaussian distribution for the B meson signal and a 
nearly flat two dimensional distribution for background. The resolutions in AE and Mes can be 
different for different decay modes. 

The intermediate state particles in this analysis are J/ifi (ee,fi/j,), <ft (K + K~), 77 (77, 7r + 7r~7r°), 
rf (77 (77) 7r + 7r~), 7r°(77),and Kg (7r + 7r~). All of the intermediate state particles are selected in 
mass windows, which are listed in Table |l]. The J/tfi —* ee decay has a slightly asymmetric mass 
window to include the radiative J/ip tail. Since B° — > J/ijjrj and B° — > J/i^r] 1 are pseudoscalar 
decays into a vector and a pseudoscalar, the distribution of the helicity angle]] of the lepton, 9 Lepton, 

The lepton helicity angle is defined as the angle measured in the J/ip rest frame between the direction of the 
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from the J/t/j is proportional to sin 2 9i epton . Hence an additional cut of |cos 9L e pton\ < 0-8 is applied 
to reject continuum and other backgrounds. For the 77 candidates, a ir° veto is applied where the 
candidate is rejected if either of the associated photons can be combined with another photon in the 
event to form a 77 mass within 20 MeV/c 2 of the 7T° mass. Also for the mode B° — ► J / ^{'yy) the 
77 candidate is rejected for asymmetric decays with cos 6!^ < 0.8, where 6^ is the photon helicity 

angle in the 77 rest frame. The rf — > 77 (77) tt + it~ candidate uses the same 77 selections, including the 
7T° veto. The mass of Kg candidates is taken at the closest distance of approach between positively 
and negatively charged tracks. 



Table 1: Mass windows used in selection of intermediate particles. 



MODE 


Mass Range (GeV/c 2 ) 




— > e + e 


2.95 < M e+e - < 3.14 


J/V 




3.06 < M(ji + fi-) < 3.14 




K + K- 


1.004 < M (K+K-) < 1.034 




-> 7T + 7T~ 


0.489 < M (it + it-) < 0.507 


7] -» 


77 


0.529 < M (77) < 0.565 


77 -> 


7r + 7T _ 7r° 


0.529 < M (vr+vr-vr ) < 0.565 


v'- 




0.938 < M (V + < 0.978 


7T°- 


♦ 77 


0.120 < Af (77) < 0.150 



An additional requirement is applied to separate and remove two-jet-like continuum events from 
more spherical B meson decays. The thrust direction of the B meson candidate and thrust direction 
of the recoiling other tracks in the event are calculated. Typically, 0t, the angle between these two 
directions is uncorrelated for BB events and peaked at cos 9t = ±1 for continuum events. The 
thrust angle requirement for the B decays is |cos#t| < 0.8. 

The PID criteria are listed in Table ^ mode by mode. The PID requirements for some modes 
are slightly more stringent for background rejection. 



Table 2: Particle identification requirements for each decay mode. 





J/tp — > e + e 


J/tp -> fi+fi 


cj) -> K+K' 3 rd K* 


j/i)(j)K + 


VTE + LE 


VLM + LM 


VTK + LK NP 


J/^K° S 


VTE + LE 


VLM + LM 


VTK + LK 




VTE + VTE 


TM + TM 


VTK + LK 


J/tpT] 


VTE + VTE 


TM + TM 




J/H 


VTE + VTE 


TM + TM 





negative charged lepton and the direction opposite to the parent B meson. 
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3.2.1 B° -> J/^cj) Mode 



This mode combines the J/ip and <p based on the selection described in the previous section. The 
resulting scatter plot of AE versus Mes is shown in Figure || (left top). The signal region is shown 
on the figure, and it is defined by 5.272 < M ES < 5.288 GeV/c 2 and -0.057 < AE < 0.057 GeV. 
The left bottom (right) plot shows the projection onto the Mes (AE) axis for events that satisfy 
the AE (Mes) requirement for the signal region. The curve overlaid on the Mes projection in 



this and the following figures is the sum of an ARGUS function [12| to model the combinatoric 
background and a Gaussian, where the Gaussian contains the background peaking in the signal 
region as well as the signal itself (see Section Q for details). Statistically there is no significant 
signal for B° — > J/if)(j). An upper limit on the branching fraction is described in the next section. 
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5.3 -0.3 
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Figure 3: AE vs Mes (left top), Mes projection in AE signal region (left bottom), and AE 
projection in Mes signal region (right) for B° — > J/ijxj). 



3.2.2 B -> J/ixj)K + and Jj^4>K% Modes 

In this mode we combine the J ftp and <fi candidates described above with a charged kaon or 
Kg candidate. The resulting scatter plot of AE versus Mes is shown in Figure |4] (left top) for 
B + — * J/ r ip(j)K + . The signal region is shown on the figure, and it is defined by 5.272 < Mes < 
5.288 GeV/c 2 and -0.057 < AE < 0.057 GeV. The left bottom (right) plot shows the projection 
onto the Mes (AE) axis for events that satisfy the AE (Mes) requirement for the signal region. 
The corresponding plots for B° — > J/ip(pKg are shown in Figure |5[ The branching fractions are 
determined in the next section. 

3.2.3 B° -> J/ipr) Mode 

For this mode, we combine a J/ip candidate with an r\ candidate in the final states 77 or 7r + 7r _ 7r . 
The resulting scatter plot of AE versus Mes is shown in Figure || (left top) for the 77 mode and in 
Figure [7| (left top) for 7r + 7r _ 7r° mode. The left bottom plot and the right plot show the projections 
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M FS (GeV/c 2 ) AE (GeV) 



Figure 4: AE vs Mes (left top), Mes projection in AE signal region (left bottom), and AE 
projection in Mes signal region (right) for B + — > J/ip(fiK + . 




M FS (GeV/c 2 ) AE (GeV) 



Figure 5: AE vs Mes (left top), Mes projection in AE signal region (left bottom), and AE 
projection in Mes signal region (right) for B° — > J/ipcfiKg. 



onto Mes an d respectively. The signal region is defined by 5.27 < M ES < 5.29 GeV/c 2 and 
-0.1 < AE < 0.1 GeV for the 77 mode, and 5.27 < M ES < 5.29 GeV/c 2 and -0.072 < AE < 0.072 
GeV for 7r + 7r~7r mode. No statistically significant signal is observed. Upper limits on the branching 
fractions for these modes are described in the next section. 
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M ES (GeV/c ) 
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Figure 6: AE vs Mes (left top), Mes projection in AE signal region (left bottom), and AE 
projection in Mes signal region (right) for B° — > J/tprj( r yy). 
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Figure 7: AE 1 vs M^s (left top), Mes projection in AE signal region (left bottom), and AE 
projection in Mes signal region (right) for B° — > J/ip7](ir + ir~7r°). 



3.2.4 B° -> J/V'??' Mode 

In this mode we combine J /ip and 77' candidates. The resulting scatter plot of AE versus Mes is 
shown in Figure I (left top) for B° -> J/VV- The signal region is defined by 5.27 < Mes < 5.29 
GeV/c 2 and -0.1 < AE < 0.1 GeV. The left bottom (right) plot shows the projected M E s (A£) 
distribution for events that satisfy the AE (Mes) requirement for the signal region. There is no 
significant evidence for B° — > J/ipr/'. An upper limit on the branching fraction is determined in 
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the next section. 




5.2 5.3 -0.3 0.3 



M ES (GeV/c 2 ) AE (GeV) 

Figure 8: AE vs Mes (left top), Mes projection in AE signal region (left bottom), and AE 
projection in Mes signal region (right) for B° — > J/iprf . 



4 Efficiencies, Backgrounds and Systematic Uncertainties 

The efficiencies for each mode are determined by Monte Carlo simulation where three-body phase 
space is assumed for the three body modes (J/ipcf)K + , J/ip<f>Kg), two-body phase space for the 
vector- vector mode (J/ip<ft), and helicity amplitude matrix elements for vector-pseudoscalar modes 
(J/tprj, J/ipr]'). The statistical error due to the number of Monte Carlo events is included as part 
of the systematic error. 

The background in the Mes distributions can be described by an ARGUS function for the 
combinatoric background, plus a Gaussian function for the peaking background. The combinatoric 
background, denoted Nargus, is due to continuum events, BB events with at least one J/ip, and 
BB events without a J ftp. The peaking background, denoted Nj/^_Q auss , comes only from BB 
events with a J/tp. The shape of the ARGUS term is determined mode by mode by fitting an 
ARGUS function to the Mes distribution from a special set of events in the data where the J/tp 
is replaced by a fake J/ip. The fake J/ip is selected with identical selection criteria in each mode 
except for logically reversing the lepton identification. This provides a large sample in each mode 
whose Mes distribution can be fitted and represents the ARGUS shape. The normalization of 
the combinatoric background for each mode is obtained from a fit to the Mes distributions in 
the AE signal region of the on-peak data. The integral of this function in the signal region is 
Nargus- This method of determining Nargus has been checked with Monte Carlo simulation, 
off-peak data and AE and J ftp mass sidebands from on-peak data. The peaking background 
N j lip-Gauss is determined from a sample of Monte Carlo BB events that is normalized to the 
equivalent data integrated luminosity and contains at least one decay of J/tp — > leptons. The 
Mes distribution from this sample is fit with an ARGUS function and a Gaussian in the AE 
signal region where the normalizations are allowed to vary. The number of events in the resulting 
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Gaussian fit is the Monte Carlo estimation of the peaking background Nji^_ Q auss . The sum of 
N ARGUS plus Nj/^_ Gauss gives n&, the total number of background candidates in the signal region 
and its error, cr&. The combinatoric background is by far the dominant background in all modes 
except the B° — > J/-0^(7r + 7r _ 7r°) mode, where the peaking component reaches ~ 20% of the total 
background. 

The following sources of systematic uncertainty are considered. 

• Uncertainty in the number of BB events (column labeled AN BB - in Table [3j). 

• Uncertainty from secondary branching fractions (column labeled SBF in Table ||) . 

• Monte Carlo statistical error (column labeled MC in Table ||). 

• Uncertainties in PID, tracking efficiency and photon detection efficiency (column labeled 
PidTrkG in Table |). 

• Variations in the event selection criteria (column labeled EvtSel in Table |3|). 

• Background parameterization (column labeled BkgdP in Table ||). 

The secondary branching fraction uncertainty combines all errors from the Particle Data Group 
(PDG) Jl^] for each mode. The fractional uncertainty in N BB is 1.6%. The uncertainty from PID, 
tracking efficiency and photon detection efficiency is based on the study of the control samples. 
The uncertainty due to event selection includes varying all event selection criteria by a reasonable 
amount and determining the effect on the branching fraction. The uncertainty from background 
parameterization is estimated by using AE sideband information. The largest systematic error 
comes from varying the event selection criteria and no single variation dominates this systematic 
in any mode. 

The total systematic error combines all these separate errors in quadrature mode by mode. The 
individual systematic uncertainties are listed in Table ||. 



Table 3: Systematic error summary. 



Mode 


^ BB 


SBF 


MC 


PidTrkG 


EvtSel 


Bk f 


;dP 


Total (<jt) 




1.C 


i% 


2.2% 


1.6% 


6.7% 


11.7% 


12.( 


)% 


18.3% 


J/^K+ 


i.e 


i% 


2.2% 


1.6% 


8.2% 


11.5% 


5.9 ( 


Vo 


15.6% 


J/i><t>K° s 


i.e 


;% 


2.2% 


2.1% 


8.3% 


14.8% 


1.9! 


Vo 


17.5% 


J/tpV (77) 


i.e 


;% 


1.8% 


1.6% 


2.9% 


14.3% 


6.9 ( 


Vo 


16.4% 


J/lp7] (37r) 


i.e 


i% 


2.4% 


2.2% 


7.7% 


13.9% 


8.0 ( 


Vo 


16.5% 




i.e 


i% 


3.8% 


4.6% 


5.7% 


11.7% 


7.1 ( 


Vo 


16.1% 



5 Branching Fractions and Upper Limits 

The branching fraction determination uses a simple subtraction of events in the signal region. The 
number of signal events is n s = iiq — rib, where the term no is the number of data events in the 
signal region, and rib is the total background described in Section 
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The modes J/vp(fiK + and J/ipcfiKg have significant signals: J /ilxf)K + is 3.1 statistical standard 
deviations from zero, while J/ip(j)Kg is 2.7 statistical standard deviations from zero. The calculated 
branching fraction is based on the Monte Carlo efficiency, n s , N BB -, and the secondary branching 
fractions for the J/ip, 4>, and Kg from PDG (T^|. The results are summarized in Table ^including 
the total summed background events in the signal region. The first error is the statistical error, and 
the second error is the systematic error <jt taken from Table [3|. The derived result for B° — > J/ipcpK 
is also shown in Table ||. 



Table 4: Branching fractions for J ' /il)4>K + ', J/t/jcf)K s and the derived result for J/ip(j)K . 



Mode 


Efficiency 


n 


n s ± a(n s ) 


n b ± a b 


Branching Fraction 


J/ijj(j)K + 

J/i>(f>K s 

J/4>(j)K 


10.6% 
8.6% 


23 
13 


15.2 ±4.8 
9.7 ±3.6 


7.8 ±0.6 
3.3 ±0.4 


( 4.4 ± lA(stat) ± 0. 7 (sy st)) xl0~ b 
( 5.1 ± 1.9(stat) ± 0.9(syst))xl0- 5 
(10.2 ± 3.8(stat) ± 1.8(syst))xl0" 5 



For the modes with no signal or weak statistical evidence (J/t/jcj), J/ipr], J /iprf) an upper limit 
is set. The upper limit method uses the number of data events counted in the signal region, no, n b , 
and its error a b (described in Section |4|), in the signal region and the total systematic uncertainty 
(IT (%) from Table § Once we obtain no, n b ± cr&, and cjt, then we assume these two uncertainties 
((7b, ax) are uncorrelated and Gaussian, the upper limit -/Vg % is obtained by folding the Poisson 
distribution with two normal distributions for these two uncertainties and integrating it to the 90% 
confidence level. We list the variables in Table |B| to obtain, iVgo%, the number of events for a 90% 
upper confidence limit. Then using the upper limit -/Vgo%, the efficiency and N B g, we determine 
the resulting upper limits on the branching fractions which are also shown in Table ||. 



Table 5: 90% upper confidence limits. 



Mode 


Efficiency 


no 


n b ± a b 


(7 T (%) 


^90% 


90% C.L. Upper Limit 


B -» J /^4) 


12.1% 


1 


0.3 ±0.2 


18.3 


3.70 


< .95 x 10" 5 


B -> J/'ipr] (77) 


15.5% 


8 


1.7 ±0.4 


16.4 


11.8 


< 3.0 x 10~ 5 


B -» J/tp7] (tt+tt-tt ) 


8.7% 


4 


1.5 ±0.9 


16.5 


6.86 


< 5.2 x 10~ 5 


B — > J/ipr] combined 












< 2.7 x 10~ 5 


B j/vv 


2.5% 





0.5 ±0.3 


16.1 


1.84 


< 6.4 x 10" 5 



6 Conclusions 

We observe evidence for B — > J/%jxj)K in two modes and determine the branching fractions 
B(B -» J/il)(i)K + )= (4.4 ± lA(stat) ± 0.7(syst)) x 10" 5 and B(B -» J /xfj(f)K s )= (5.1 ± 1.9(stat) ± 
0.9(syst)) x 10~ 5 . The branching fraction for B — > J/ip(f)K is consistent with CLEO results Q. 
Upper limits have been determined for the modes B — > J/ip(fi, J/tpTj, and J/iprj'. However, the 
two B — > J/tpr] upper limits in Table || would correspond to a combined branching fraction of 
(1.6 ± 0.6(stat.) ± 0.2(syst.)) x 10 -5 , which is comparable to the B — > J/ipir branching fraction. 
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